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We consider inelastic collisions between atoms of different kinds as a 
potential source of photons in the sonoluminescence phenomena. We esti- 
mate the total energy emitted in one flash and the shape of the spectrum 
and find a rough agreement between the results of our calculation and 
the experimental data. We conclude that the atomic collisions might be 
a candidate for the light-emitting mechanism for sonoluminescence and 
discuss the implications. 
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The sonoluminescence is a mechanism of the conversion of sound energy into the 
energy of light emitted in a picosecond flash. This phenomenon was discovered a long 
time ago (the emission of light from the water was observed in but only recently has 
attracted much attention because it becomes possible to trap a single bubble of a gas 
mixture in a sound field (this is so called Single Bubble Sonoluminescence — SBSL 
or just SL). The detailed review of the topic can be found in ||. The mechanism of 
sonoluminescence has not been definitely understood so far. There are several possible 
explanations, each of them possessing its own difficulties. Let us briefly review some 
of the models. 

One of the most promising mechanism is the Bremsstrahlung radiation from ionised 
regions [|]|| . These regions are created by shock waves formed during the collapse of the 
bubble. This model succeeds, in particular, in prediction of the total radiation of energy 
per flash. However, the presence of the plasma would imply the sensitivity of this light- 
emitting mechanism to the external magnetic field |J , which is not observed. Another 
possible explanation is given by the model of blackbody radiation J7|]. The spectrum 
of the SL is well fitted by the blackbody radiation spectrum of temperature 25000 K 
(in water of temperature 22°C M). However, the time dependence of SL spectrum is 
independent of the color of the emitted light, which contradicts the blackbody radiation 
model expectations 0. Other models refer to purely quantum effects. These are: the 
"dynamical Casimir effect" suggested by Schwinger [|1(J and the Unruh effect elaborated 
in more details by Eberlein |Tl| . Unfortunately, the dynamical Casimir effect gives 



ambiguous predictions depending on the renormalization procedure that one chooses, 
while the Unruh effect gives the power of radiation too low by some orders of magnitude. 
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Yet another model describes the SL radiation as a radiation emitted from molecular 
collisions jEJ. This model properly predicts the shape of the spectrum, however it tends 



to overestimate the magnitude of the emitted energy if one uses the realistic values of 
the radiating volume and the emission time. Besides, recent developments suggest that 



due to sonochemical reactions mainly noble gases are left inside the bubble [13|, which 
if true, may be troublesome for the mechanism of SL. Since none of the mechanisms is 
fully satisfactory we may, in search for an alternative candidate, address the question 
if the inelastic collisions between atoms may be the source of light in sonoluminescence. 

Thus we consider in this article the exclusive reaction A\ + A2 — > A\ + A2 + 7 
of photon production in which the atoms remain in their ground state. There are 
interesting features concerning these reactions which encourage to study them in more 
detail: 

(i) The featureless spectrum of SL suggests that the emission of light from excited 
atoms do not play a crucial role in the phenomenon. 

(ii) The spectrum of photons emitted in atomic collisions can behave similarly to that 
of SL. Indeed, the emission of photons with a wavelength greater than the size of the 
atom is suppressed due to the mutual cancellations between the contributions of the 
nucleus and the electronic cloud (c.f. the Rayleigh scattering process). The short 
wavelength photons are absent because of the cut-off on the energy available in the 
collision. 

In the case of SL from an air bubble in water the analysis of the scattering of 
the oxygen or nitrogen atoms on the atoms of argon should be performed. The exact 
calculation for this process is extremely involved thus, using a simplified model, we 
estimate the order of magnitude of the cross-section and the shape of the spectrum. 
Hence, let us start from the process of scattering of a hydrogen atom in the Coulomb 
field of an infinitely heavy source, accompanied by emission of light. After obtaining 
the cross-section we shall argue that this elementary process has similar features to the 
atomic collision and we shall use the results to estimate the number of photons in the 
SL flash. 

We employ the standard non-relativistic hamiltonian whose leading part H con- 
tains the operators of the kinetic energy and the electron-proton interaction. The 
perturbation Hi in the Coulomb gauge, after neglecting A 2 terms, take the following 
form: 

gi(r..r,) = - eA(T ' )P * + eA(T ' )P * - + gg, (1) 

m e m p r e r p 

where r e (r p ) and p e (p ) are the position and momentum of the electron (proton), the 
source of the potential is situated in the origin and carries the electric charge of Ze, 
Ak(r) = Sk^y=a^ k is the vector potential of the emitted plain wave with momentum k 
and the polarisation vector given by The matrix element of the creation operator 
in the box is given by (7(5, fc)|a+ fc |0) = y^-- The eigenstates of H Q are given in the 

position space by the wave functions ^p >n (R, r) = 1/y/V exp(— iPR) ip n {r) with P 
being the atom momentum, ip n (r) — the electron-proton wave function (for a bound 
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or continuum state n) in the centre-of-mass frame and r = r e — r p , R = ar p + j3r e , 
a = m p /(m e + m p ), (3 = m e /{m e + m p ). For future use, we also introduce the mass 
of the atom Mh = m e + m p . We shall calculate the amplitude Ai of the transition: 
*&p,is —> ^P',isl(k) in the second order of perturbative expansion which gives the first 
non-vanishing contribution. In this approximation the amplitude reads 

M = J- (^M^( fc )l^il m )( m l^il^.ig) ( 2 ) 
m E P:1S — E m + ie 

where m runs over all the possible intermediate states. There are two groups of the vir- 
tual states depending on the photon content. The first one contains only the electron- 
proton states moving with the total momentum P" , which may be either virtual bound 
states or virtual continuum states. We shall denote the members of this group by 
\^fp" } n). The other group includes the states \^p» t n', 7(&)) built of a photon of energy 
k and momentum k accompanying states |\I/p». n ) in the electron-proton sector. It is 
straightforward to observe, that the sum (|2]) may be expanded into two parts corre- 
sponding to final (FSR) and initial state (ISR) photon radiation amplitudes. ^From 
the kinematics of the process and the conservation of momentum it is easy to find 
that the energy denominators take the following form: -Dfsr — Eig — E n + k and 
.Disr — Eis — E n — k where, making the approximations, we have employed the hier- 
archy k < P ~ P' < M H - 

We have checked, that the contributions to the amplitude Ai of the intermediate 
states in which the electron and proton form a IS state are negligible, due to a^k 2 
suppression of the photon emission matrix elements. On the other hand, the spectrum 
of hydrogen atom has a large gap between the ground state and the first excited state 
and the excited bound states lie close to each other. It also may be proven that the 
intermediate continuum states with energies E n > \Ei$\ give a sub leading contribution 
to the amplitude. It seems therefore quite natural to approximate the energy difference 
Eis — E n by a characteristic energy, say E\s — E n ~ — AE = —\Eis\. Of course, this 
approximation works only for k significantly smaller than the AE. Now, we can employ 
the completness relation for the ^p tn states in the electron-proton sector and perform 
the sum over P, n to obtain [] 

M ^ ^AE—k + ^AE~+k ' (3) 

where H$ and He are the components of Hi responsible for the atom Coulomb scat- 
tering and for the photon emission respectively. After performing the necessary inte- 
grations and taking into account the relation Sk k = we get the following expression 
for the amplitude: 



This approach is in the spirit of the classical approximation applied in the calculation of 



van der Waals potential [14] 
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{ AE 2 - k 2 

2AE 
AE 2 - k 2 



/ 2n AnZe3 
~ V kV 3 (P' + fc-P) 2 * 



P-f{(3q + k) + ^-f{(3q) + —f{aq) + —f(aq - k) 

m e m p m e m p 



+ 



Sk(P + P') lf(aq) - f(aq + k) + f(J3q) - f(J3q + fc)]} (4) 



where q = P' — P and f(l) = (1 + a^l 2 /4)~ 2 is the electric form-factor of the ground 
state of the hydrogen atom of the Bohr radius given by a . As expected, there occur 
cancellations between the form-factors, and the amplitude vanishes for vanishing k. We 
also notice that the terms proportional to £k{P + P') are suppressed by a small factor 
aok as compared to term containing e^q thus we retain only the later. Furthermore 
we may drop the l/m p terms. 

A more subtle problem arises when comparing the relative importance of the con- 
tributions ^f(/3q + k) and -^-f(aq). The answer seems to depend on the details of 
the process in the physically relevant parameter space. For example, if the energy of 
the collision falls between 5 and 10 eV then the first term dominates for k > 3 eV while 
for k < 3 eV the contribution of the other term is more important giving rise to an 
additional "red" peak of spectral density. However, if we take an atom with the mass 
of 16 Mh and the Bohr radius ao (modelling the oxygen atom) this peak would move to 
energies k < 0.5 eV i.e. to the experimentally inaccesible region. Besides, the ^-f(aq) 
term is very sensitive to the details of the charge form-factor of the scattering atom. 
On the other hand, the —f((3q + k) part is universal i.e very weakly dependent on 
the details of the atom structure. In what follows, we shall focus only on this universal 
contribution of the amplitude which is relevant for the shape of the spectrum in its 
large k part. By this neglection we underestimate the number of the emitted photons. 

Thus we get the estimate 



= / 2tt AnZe 3 k e k g 

V kV 3 M H AE 2 - k 2 q 2 [ ' 

/,From this amplitude, after performing the standard integrations over the phase space, 
one gets the differential cross-section for the photon emission from a collision at the 
energy E: 

da _ 8Z 2 e 6 k 3 / P + P \ 

dk ~ 3M H E (AE 2 - k 2 ) 2 ° g \P~^P ) ' U 



where P = ^2M H E, P = y2M H (E — k). For the purpose of the numerical calcula- 
tions, we take AE = 13.6 eV. 
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Let us point out, that the obtained cross-section corresponds to a subprocess in 
which the nucleus scatters off the Coulomb field and the electron cloud radiates a pho- 
ton when recombining around the scattered nucleus. The other contributions to the 
amplitude were shown to be subleading and neglected. The spectral density resulting 
from this cross-section rises as k 4 for small k i.e faster than this observed in the SL 
experiments. For k closer to the kinematical cut-off, the initial k 4 dependece is mod- 
iffied and one gets approximately quadratic behavior over wide range of k which is still 
steeper that it follows from the experimental data. However, let us remind, that in our 
picture, the SL spectrum would be given by a convolution of the spectra of individual 
collisions with the (unknown) distribution of the collision energy, so the two spectra 
may differ from each other. 

After integration of the differential cross-section (^) over the energy, we get the 
total cross-section <Jtot{E) growing with the collision energy E approximately as E v 
with v ~ 3.5. For the choice of parameters relevant for the hydrogen atom, Z = 1 and 
a typical collision energy of 7 eV (the choice is suggested by the observed cut-off on 
photon spectrum) we obtain cr = a tot {7 eV) = 1.3 • 1CT 31 m 2 . Taking this number and 
the parameters implied by the known facts concerning the bubble dynamics we may, 
after some modifications, make a crude estimate of the number Nf of photons produced 
in one flash. We assume, that the bubble is filled with a mixture of atoms of a noble gas 
and atoms of other element (oxygen or nitrogen) | with atomic numbers Zn and Zo, 
masses and Mq and concentrations and no correspondingly. The predictions 
for collisions of an atom of noble gas with an atom of, say, oxygen may be formulated 
on the basis of the previous calculation. Namely, for the momentum transfers relevant 
for SL the condition a^q 2 ^> 1 is fulfilled and the charge form-factors for the electrons 
suppress their contribution to the scattering amplitude — the electrons cannot absorb 
momenta much larger than their average momentum in the atom. Then the nucleus of 
the noble atom acts as a bare source of the Coulomb field. Furthermore, the amplitude 
of the electromagnetic radiation from the noble atom is smaller then from the oxygen. 
It is caused by the smaller polarizability of the former atom which is governed by the 
relevant ionisation energy. Thus the cancellations which occur due to the destructive 
interferece do not reduce the cross-section substantially in contrast with the case of a 
collision of two objects of the same kind. This feature was reflected in our model by 
neglecting the radiation from the source of the Coulomb field. Of course, to complete 
the analogy, we should use in the formulae the reduced mass of the two atoms instead 
the mass of the hydrogen atom and include the modification of charges by Zm for the 
source of the field and Zo for the scattering atom. The charge Zo enters in fourth power 
since it contribute both to scattering and the radiation. Therefore we approximate Nf 
by the following expression: 

Nf^—vr MqMn Z n Z a n N n , (7) 



2 At the temperatures and densities predicted by a shock-wave model the diatomic molecules 
are dissociated. 



5 



where R s is the radius of the hot gas region, r = 100 ps is the light emission time 
and v ~ 10 4 m/s denotes the relative velocity of colliding atoms. Let us focus 
on the argon (Zjy = 18)-oxygen (Zo = 8) process, which may be relevant for the 
SL of air bubble in water. As a first guess we take for R s the minimal radius of the 
bubble i.e. about 0.5 /im and assume for the concentrations no + un — 600no in 
accordance with the measured compression factor @, with n = 2.7 • 10 25 m" 3 being 
the concentration of the ideal gas in the normal conditions. We obtain Nf ~ 5 • 10 5 
for no = n^-, in a reasonable agreement with the experiment. In the alternative shock- 
wave description of gas dynamics lower values of the radius are prefered, however the 
increase in concentration caused by the additional compression may easily compensate 
the decrease of the reaction volume. In fact, the true concentration in the radiating 
region is not measured, and the models provide us only with ambigous predictions. 
The limitting value of the concentration is ~ l/(8ag) = 10 30 m~ 3 thus there is still 
some room here and the constraints imposed on the cross-section (|7|) by the experiment 
are not very stringent. It is also very probable, that the proposed mechanism is not 
responsible for the emission of all photons but rather supplements the list of processes 
studied in ||. In this case, including it could, in particular, improve the shape of the 
spectrum obtained in H. 



We can also take into account that, as follows from [13], the SL bubble is filled 
mainly with a noble gas, and may contain only some admixture of oxygen (which may 
be continously provided from water) by modifying the ratio no : n^, keeping + % 
fixed. Thus for 10% of the oxygen in the gas mixture the number of photons in one 
flash drops to about 2 • 10 5 . 

In conclusion, we propose a novel light-emitting mechanism in sonoluminescence, 
in which the photons are radiated from atoms disturbed by collisions. The number 
of produced photons has the proper order of magnitude, however the resulting photon 
spectrum comes out somewhat to steep. The spectral density is featureless and univer- 
sal. Our model of SL requires the gas temperature to be approximately 30 000 K, and 
the concentration of the order of 500 — 1000 no- The bubble should contain atoms of 
two different g noble gas and a gas with smaller ionisation energy e.g. oxygen, 

however one of the gases may appear at much smaller concentration than the other. 
The emitting process is insensitive to external magnetic field and does not require the 
presence of molecules which are expected to dissociate in the gas temperature sug- 
gested by the photon spectrum. Although the model we have elaborated is based on 
the simplified assumptions we expect that the obtained predictions estimate correctly 
the order of magnitude of the cross-section. Therefore this gives the strong motiva- 
tion to further study of the subject which would probably require the use of numerical 
analysis if one intends to cover all the details of the process. 
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